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1. INTRODUCTION 


The vast majority of composite laminates in use today are fabri- 
cated so that the material properties are symmetric with respect to the 
laminate's geometric midsurface. To have symmetric material properties 
a laminate must be fabricated so that for every lamina above the geometric 
midsurface with specific material properties, thickness, and fiber 
orientation, there is a lamina below the midsurface with identical 
material properties, thickness, and fiber orientation. A laminate so 
constructed will not exhibit bending-stretching coupling in its elastic 
response to external loads. Coping with analyses which must include 
bending-stretching coupling effects can be difficult. Thus to simplify 
the analysis of composites, laminates a^e fabricated with midsurface 
symmetry of material properties. However, to make a laminate symmetric 
only to simplify the analysis is a penalty on the performance of the 
composite. In practice, symmetry is generally achieved by adding a 
sufficient number of layers to the laminate. From a strength or stiff- 
ness measure, these lamina may not be necessary. However, adcing the 
lamina most certainly increases the weight of a structure. It may be 
that unsymmetric laminates are actually desirable. Experimental evidence 
has shown that unsymmetrical ly laminated composites, which are flat at 
their elevated cure temperature often cure to a shape closely approxi- 
mated by a right circular cylinder. The out-of-plane warping of the 
cured laminate is a direct result of bending-stretching coupling. 

However, the cylindrical characteristic is contrary to the predictions of 
classical lamination theory LI, 2]. The classical theory predicts the 
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cured shape of all unsymmetrically laminated composites to be saddle- 
like, having their principle curvatures of opposite sign but not neces- 
sarily of the same magnitude. For some urisymmetric laminates, however, 
the saddle-shape does appear. It was felt by many investigators that 
the observed cylindrical shapes were due to problems with curing, 
moisture absorption, or other unwanted effects in the fabrication process. 
However, a systematic experimental study [3] has shown that the phenomena 
is repeatable and whether a cured unsymmetric laminate assumes a cylin- 
drical shape or a saddle shape depends on the thickness of the laminate 
compared with a characteristic in-plane length. Thin laminates cure 
into the shape of a right circular cylinder while thicker laminates 
conform, at least qualitatively, to the saddle-shape predictions of the 
classical theory. This characteristic of thin unsymmetric laminates is 
felt to be very important because of the large number of components on 
aircraft which are cylindrical in nature. In addition, fabrication of 
these cylindrical shapes is not much more difficult than the fabrication 
of flat laminates. No expensive mandrels or forms are needed; the 
laminate simply cools to the cylindrical shape. 

Recently a theory has been developed which accounts for this behav- 
ior of unsymmetric laminates [4]. The theory is a nonlinear extension 
of classical lamination theory and at present is limited to the analysis 
of laminates having their lamina oriented at either 0° or 90° relative 
to some global coordinate system. However, the predictions of the 
theory agree well with the limited amount of experimental data avail- 
able. The theory indicates that the radius of curvature of the cylinder 
is a function of the lamina orientations (i.e. [90^/0^]y, [90^/0^ ]-p > 
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etc.) and the thermo-elastic properties (E-j, Eg, ot-j, otg* •••) the 
individual lamina. The theory is being extended to account for arbi- 
trary lamina orientations. It appears, even at this early stage of the 
development of the theory, that this phenomena should be strongly consid- 
ered as a way to fabricate cylindrical components from composite mate- 
rials. The next question is: How useful are these cylinders? Do they 

buckle easily? Do they have unusual vibration characteristics? How do 
they respond to point loads, pressure loads, etc.? Perhaps any unusual 
response characteristics can be off-set by the ease of fabrication 
considerations. To answer the just-mentioned questions, a rational 
step-by-step series of analyses of un symmetrically laminated cylinders 
is in order. These analyses should range from linear analyses to nonlin- 
ear analyses, static and dynamic, steady-state dynamic and transient 
dynamic, prebuckling and postbuckling. The analyses should precede in 
such a way that potential problems are first examined at a simple level, 
say, with linear analyses, and then if no problems appear, cylinders 
should be examined at a more complicated level, say, wi' n nonlinear 
analyses. At all times the behavior of the unsymmetric laminates should 
be compared with a symmetric counterpart, its competitor based on current 
thinking. This paper reports on the first step in examining the utility 
of cylinders fabricated from unsymmetric lay-ups. The paper examines 
the static structural response of an unsymmetrically laminated long 
cylindrical panel subject to a line load acting inward. A line load, 
though simple, is felt to represent one of the more severe types of 
static loadings a laminate can encounter. The analysis presented is a 
geometrically linear and elastically linear analysis. 
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The report begins in the next section by describing the geometry 
and the nomenclature of the cylindrical panel. The section continues by 
describing the kinematics of the assumed deformation. From the kine- 
matics such fundamental quantities as strain, cross-section rotation, 
and curvature are established. Next, the nomenclature associated with 
the applied load is discussed and the internal force resultants are 
defined. Newton's second law is applied to a small section of the panel 
to establish static equilibrium equations in terms of the internal force 
resultants and the applied load. Since the panels studied nave only 0° 
and 90° lamina, the constitutive behavior of the cylindrical panel is 
represented by an orthotropic stress-strain relation. The out-of-plane 
stresses are assumed to be negligible and so the constitutive behavior 
ultimately reduces to the familiar A, B, and D matrices. In the study 
of these unsymmetrical ly laminated cylindrical panels it is shown to be 
convenient to define a reference surface other than the laminate's 
geometric mid-plane. The location of the reference surface through-the- 
thickness is related to the material properties of the panel. This 
aspect is discussed in detail. After the discussion of the constitutive 
behavior, the boundary conditions and transitions conditions at the 
point load are examined. This next section concludes by establishing 
the differential equations which govern the deformation behavior of the 
panel . 

The third section of this report is devoted to solving the governing 
equations. The boundary and transitions conditions are applied and 
closed-formed expressions for the deformation of the panel, as a function 
of spatial position along the panel, are obtained. The forrth section 
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of the paper presents numerical results regarding deflections, stiff- 
nesses, and strains of several representative panels. The numerical 
results for unsymmetrically laminated panels are examined in the context 
of results for a symmetrically laminated panel which use the same amount 
of material. This type of comparison is done to illustrate any advan- 
tages or disadvantages of unsymmetrically laminated panels. The formula- 
tion and results presented herein are valid for arbitrarily deep cylindrical 
panels. A formulation and numerical results are presented for shallow 
panels. Comparisons between the deep theory and ths shallow theory are 
made. The main purpose of this comparison is to establish whether 
deepness of cylindrical composite panels is a purely geometric property, 
as it is for metal panels, or whether it involves material properties as 
well. This latter idea has been suggested by some investigators. 

Finally, the fifth section reports on some experimental results. 

The force deflection characteristics of a [90^/0^J T and a [ (0/90)2] s 
T300/5208 graphite-epoxy cylindrical panel are compared with theoretical 
predictions. The purpose of this comparison is to determine if predicted 
differences between symmetric and unsymmetric panels are observable. 

The experimental set-up is briefly described as is the correlation 
between the theory and experiment. 


4 


2. PROBLEM DEFINITION 


Geometry 

Figure 1 illustrates the basic notation associated with the cylin- 
der geometry. The undeformed circular cylinder has a center of curvature 
at 0 and has a semi -opening angle 0 . The laminate has thickness h. 

There are two surfaces within the laminate to be considered. Normally, 
analyses of laminated plates and shells deal with the geometric mid- 
surface of the laminates. All geometry, boundary conditions, and other 
information is referenced to this surface. To date, analyses of composite 
panels have dealt primarily with symmetric laminates. In these cases it 
is natural to consider the geometric midsurface as the reference surface 
from which to describe the panel's geometry and deformation. As will be 
shown later, duo to the elastic bending-stretching coupling inherent in 
unsymmetric laminates, it advantageous to consider another surface as 
the reference surface. Thus the analysis here will discuss both the 
geometric midsurface and the reference surface. These two different 
surfaces are illustrated in Fig. 1. The laminate midsurface and reference 
surface are a distance ' d ' apart. The quantity d, to be defined later, 
is a material property. The radius of curvature of the laminate mid- 
surface is 'a' while the radius of curvature of the reference surface is 
'R. 1 These quantitites are related by 

R = a + d. (1) 

The arch rise, denoted as 'H,' is measured at midspan and is the ver- 
tical distance from the midsurface to the horizontal line joining the 
end-points of the midsurface arc. The variable 9 is the independent 
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spatial variable. It describes the angular position along the arch, and 
it is measured positive clockwise from the midspan vertical. Position 
through the thickness of the laminate is denoted as ‘z,‘ with z = 0 at 
the laminate reference surface. 

The radial displacement of the reference surfc'e is denoted as 'w,' 
positive w being towards the center of curvature. The circumferential 
displacement cf the reference surface is denoted as ' u , ' with positive u 
being to the right. In this analysis, w = w(e) and u = u (0 ) . 

Illustrated in Fig. 1 are pinned supports, the pins acting at the 
ends of the geometric midsurface. The supports are constrained to 
remain a fixed horizontal distance apart, this distance being denoted as 
'L.' The results presented here are restricted to these pinned conditions 
but sufficient information is given in deriving the governing equations 
to allow other support conditions, e.g., clamped, to be studied. 

Assumed Deformation of the Reference Surface 

Figure 2 depicts the kinematics of the assumed deformation. Due to 
the applied load, each point on the panel displaces an amount A and the 
cross-section at each point rotates through an angle ft. In the notation 
used here, the overbar denotes a vector quantity. The unit vectors 

A A 

i and i Q denote unit vectors normal to and tangent to the undeformed 

/N A 

reference surface. Likewise, the unit vectors n and t denote unit 
vectors normal to and tangent to the deformed reference surface. The 
orientation of these vectors are all a function of position alon^ the 
panel. The location of a point on the undeformed reference surface is 
given by 


( 2 ) 
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The location of that point after deformation is located by 

r = R + A . (3) 

The displacement of the point is represented by 

A = -w(e)i r + u ( 0 ) i Q . (4) 

Using Eqs. 2-4, the extensional deformation of the reference 
surface can be determined. If dS represents the magnitude of an element 
of arc length of the undeformed reference surface and ds represents the 
magnitude of the arc length of that element when it is deformed, then by 
definition the extensional strain of the reference surface is 


- ds - dS 
0 " dS 


(5) 


Arc length in the undeformed geometry is 

dR = (Ri r (e))de = Ri 0 (e)de, (6) 

where use is made of the fact that 


dl (e) 

-4-=i e (e). (7) 

The magnitude of the arc length in the undeformed geometry is 


dS 
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dR*dR 


Rde . 


( 8 ) 


In the deformed geometry arc length is defined similarly, namely, 
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dr = 


dR dA 
de de 


de 


Ri e (e) + M (-w(e)i r (e) + u(e)i 0 (e)) 


de. 


Using Eq. 8, the notation ( )• = and the relation 


d i 0 ( 6 ) 


de r 
Eq. 10 becomes 


= -i‘ (e) , 


dr = Ri fi - w'i r - wi Q + u'i 0 - ui r de 


*K S 


e 

, \ 
w' + u 


i +( ] + ULl_W M 
r \ R / 'e 


Rd6 . 


By defining 


r r=-(^H r 0=f^t- W 


the deformed arc length can be written as 


dr = R 1 + r e) ? 8 + r r ? r] dS • 


The magnitude of the deformed arc length now becomes 

1/2 


- 0 o 1/2 

ds = -ydr*dr = [(1 + r 0 r + rj; ] dS . 


(9) 


( 10 ) 


(ID 


(12) 

(13) 


(14), (15) 


(16) 


(17) 


The extensional strain of a line element is then computed directly to be 


E e * " V (1 ♦ r e> 2 + > ' 1 • 


(IB) 




J 
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Rewriting this as 
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e« = ^1+ 2T n + r? + r 2 


'0 


■0 0 


- 1 , 


(19) 


and using the binomial expansion of the first term, the extensional 
strain can be written as 


e 0 = r 0 + \ T l + 0(r3) • ( 2 °) 

Assuming the deformations of the panel are small, 

0 < | r 0 | « 1 and 0 < |r r | « 1 , (21), (22) 


the higher order terms can be dropped. The strain-displacement relation- 
ship then becomes linear, namely, 


e 0 




(23) 


The rotation, n, and the curvature, k, at points of the reference 
surface can be related to u and w by further analysis of the kinematics. 
Again referring to Fig. 2, the tangent and normal to the deformed 
reference surface are related by the rotation to the tangent and normal 
of the undeformed surface. These relations are 

A /\ /N 

t = - sin(fi)i r + cos(n)i Q (24) 

A /\ /\ 

n = cos(fi)i r + sin(f2)i 0 . (25) 


By definition 


- f j f <, 

IPC v ■ v •• ■■ 
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2 _ dr 
* = ds * 


( 26 ) 


Using previous results, 


t = 


... dr _ 1 dr _ 


ds 1 + e« dS 
u 


t 1 + Ve + Vi 


1 + e. 


(27) 


Comparing Eq. 24 and Eq. 27, 


1 + r. 


cos(ft) = j 


+ e ft * 


sin(fi) = 


-r\ 


i + e 9 ■ 


(28), (29) 


Thus the tangent of che rotation is given by 

-r 

tan(fi) = -j “ — . 


(30) 


If the deformations are small, then in addition to Eqs. 21 and 22 being 
true, 


1 + e Q : 1 . 


(31) 


For these conditions the rotations are approximated by 


a 8 -r. = 


w 1 + u 


r R 


(32) 


The curvature of the reference surface is defined by 


dt _ ' 

Hs 8 " cn 


and 


ds = - Kt 


(33), (34) 


Again, using previous results, Eq. 33 can be expanded as 
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dt _ dt dS _ 1_ dt 

ds dS ds 1 + Eq dS 


14 o rigim&. r;: ' 

a 0F P0D.1 C >- •' 

J dt 

+“£^T de • 


Using Eq. 24 in Eq. 35 leads to 


| - cos(ft) ft' i r - sin(ft) i Q 


- cos (ft) i 


- sin(ft) ft'i. 


Regrouping results in 


dt _ 1 + a 1 ) 


cos(ft)i r - s 


in(ft)i Q i . 

9 I 


Substituting Eq. 25 in Eq. 37 leads directly to an equation of the form 
of Eq. 33, namely 


dt _ (1 + ft') : 

ds " " R (TT^T n 


Therefore, curvature is defined to be 


1 + ft' 


For small deformations 


K = 1+ £Ln K + fil 

K R R K o R ’ 


k q being the initial curvature of the reference surface. Substituting 
for the definition of ft, Eq. 32, the curvature can be related to u and 


w by 
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K a 


w" + u‘ 
' R 2 
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Deformations off the Reference Surface 

A similiar kinematic analysis can be developed for lines parallel 
to the reference surface. Figure 3 shows the deformation of a point 
which is a distance z away from the reference surface. For purposes of 
illustration, the thickness of the laminate has been exaggerated in Fig. 
3. Due to the applied load, the point displaces by an amount 3 _. The 
subscript z is used to denote off-reference-surface position. It is 
assumed that deformations normal to the reference surface, i.e. through- 
the-thickness deformations, are negligible in this analysis and that 
there are no transverse shearing deformations. Thus, in the deformed 
state, the point is still a distance z from the reference surface and 
the distance z is measured along the normal to the deformed reference 
surface. 

The extensional strain on a parallel line is also determined by 
examining arc lengths. The undeformed location of the point off the 
reference surface is 

R z = (R + z) i r . (42) 

From this, 

dR z = (R + z)dei e =^1 + f^dSi ? e . (43) 

and 

ds z = ^dR z • dR z = (1 + ZK Q )dS . (44) 


In the deformed geometry, the point is located by 
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Fig. 3 


Kinematics of Deformation 
Off the Reference Surface 
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r; 



r - r + zn . 


( 45 ) 


The deformed arc length can be determined by 


) 


dr z = (3? +z i? ) ds = (1 +ZK) t ds 


dr = (1 + ej (1 + zk) t dS. 
z u 


(46) 

(47) 


In going from Eq. 46 to Eq. 47, Eqs. (26), (34) and (5) were used. The 
magnitude of the deformed arc length is 


ds. 


= ^dr • dr = (1 + e 0 )(l + zx)dS . 


(48) 


Substituting Eqs. 44 and 48 into the definition of extensional strain, 


ds_ - d$ 
t 0 dS 


(49) 


leads directly to 

Eg = (1 + Eg) (1 + zk) (1 + ZK 0 ) _1 - 1 . 


(50) 


Expanding the first two terms and using the bionomial expansion of the 
third terms results in 


Eg = Eg + Z(1 + £g)(K - K q ) + 0(Z ) . 


(51) 


For small deformations, the extensional strains on parallel surfaces 
become, using Eq. 31, 


E 0 = E 0 + Z(K - K 0 ) ■ 


(52) 
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Parallel surface displacements can be computed in terms of reference 
surface displacements and rotations by considering the geometry of Fig. 

3. From the figure 

zi f + A z = I + zn . (53) 

Using Eqs. 4 and 25 and grouping terms produces 

I 2 = [-w + z(cos(ft) - 1 ) ] i r + [u + z sin(ft)]i 0 . (54) 


If A z is written as 

\ = -w(6,z)i r + u(0,z)i0 , (55) 

then by association in Eqs. 54 and 55, 

w(0,z) = w(6) + z(l - cos ( fi( 6 ) ) ) (56) 

u(0,z) = u(0) + z sin(C(e)) . (57) 

For small rotations, cos(fi) ~ 1, 

w(3,z) = w(8) , (58) 

and 

u(6,z) = u ( 6 ) + zft(e) . (59) 


Applied Loads, Internal Force Resultants, and Static Equilibrium Equations 
The load analyzed in this study is a line load of intensity P, with 


P in units of force per unit length along the generator direction. 
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Referring to Fig. 4, the force P is positive when directed toward the 
center of curvature and the load acts at some arbitrary location on the 
span denoted by 0 = The load is inclined from the radial direction 
by an angle a. 

The internal force resultants are illustrated in Fig. 4. The in- 
plane (membrane) force, N, is assumed positive in tension. The positive 
sense of the shear resultant, V, and moment resultant, M, are indicated 
in the figure. These resultants are assumed to act at the reference 
surface and, as in a shell theory, are defined per unit length along the 
generator or x-direction. 

To determine the governing equilibrium equations, a small portion 
of laminate is isolated. The pertinent forces acting in this portion of 
the laminate are illustrated in Fig. 5. Since this is a linear theory, 
the forces are illustrated in and Newton's Law will be applied in the 
urideformed geometry. Even though the interest here is in a line load, a 
distributed load q(6) is shown in the figure for generality. It is 
assumed spatially uniform in the x-coordinate and has units of force per 
unit reference surface area. 

A 

Referring to Fig. 5 and summing forces in the normal, i r , direction 
leads to 

- N = Rq cos(a) . (60) 

A 

Summing forces in the tangential, i 0 , direction results in 


+ v = R q sin (a) . 


(61) 


20 


original mm sy 

OF POOR QUALITY 



Fig. 4 


Applied Forces 


and Internal Force 


Resultants 
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Finally, moment equilibrium of the element requires 

+ RV = 0 . (62) 

Equations 60 - 62 are the static equilibrium equations which govern the 
behavior of the cylinder under load. 

Constitutive Behavior and the Determination of the 
Reference Surface Location 

The constitutive equations for a lamina are 

a x = ^ll E x + ^12 E 9 
°0 r ^12 E x + ^22 E 0 
T x0 = Q 66 y x0 * 


(63) 

(64) 

(65) 


where E , E Q , y xQ are the strains and a x> a Q , t xQ are the stresses. 
Because the laminates in this study are made from lamina with fibers 
either in the 0° or 90° direction relative to the x-di recti on. 


Q^g ~ i ” 1 >2. 

i 

For a panel which is long in the x-direction, it is assumed that 


E x = Y x0 = 0 * 


( 66 ) 


such that 


a x = Q 12 E 8 ’ 


(67) 


a 9 = Q 22 E 0 * 


T “0 
xy 


( 68 ) 

( 69 ) 
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For a beam it is assumed that the transverse stresses are zero, i.e. 
a = 0 in this case. Hence 

A 


Q?2 

o 0 = ^22 ^ " W g- ^ E 9 * 

y lP22 


(70) 


In either formulation the constitutive equation can be written as 


°0 = ° 22 E 6 1 


(71) 


where 


/ Q 22 for a long panel 


Q* 


22 


0?2 


Q22 ^ " Z~Z — ) f° r a bea ^- 


Q n Q 22 


It appears the factor 0] 2^1 1^22 is sma11 W1th respect to 1 for the 
practical case. For example, consider a 0° ply of graphite-epoxy. From 
Jones [1 , pp. 46-51], 


Q 11 = Q 11 = ~ 


•1 


v 12 v 21 


(72) 


Q 22 = ^22 


1 " v 12 v 21 


(73) 


^12 = Q 12 


v 12 E 2 


1 ‘ v 12 v 21 


(74) 


Thus 
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12 


On'tl 


vf 2 E 2 
T 


(0.21 )*1 .7 nisi 

SlTiiSl — 


» 0.0037485 


1T22 


(75) 


which Is, In fact, small with respect to unity. Therefore Qgg « ^ Is 
exact for a long panel and Is a good approximation for a beam or arch. 
The definitions of the stress resultants for the reference surface 


are 


N 



> dz 
+ d) 



M ° 


-I 


2a e 

4* 


dz 

d) 


(76), (77) 


Using the constltive relation, Eg. 71, the definition of parallel 
surface strain, Eq. 52, and integrating, Eqs. 76 and 77 become 

N ° ^22^8 + B 22^' " K 'o ) 


M 


^22^8 


" D 22^' 


K o! 


(78) 

(79) 


where 


(Agg.Bgg.Dgg) 


Q 22 (l ,z,z 2 ) dz 


-(jr + d) 


(80) 


For the middle surface, making the change of variable z « i - d and 


substituting into Eq. 80 yields 
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(81) 

(82) 

(83) 


in which 





h 

2 * p 

0^2 (1* z * *) dz 


(84) 


-h 

2 

are the stiffnesses for the middle surface. To uncouple the constitutive 
equations requires 


B 


22 = 0 = B 22 - 



(85) 


so 


B 22 
d = -££ 

A 22 


and 


( 86 ) 




(87) 


The barred stiffnesses are the stiffnesses typically associated with 
classical lamination theory [1, pp. 154-155]. Using the uncoupled 
reference surface, Eqs. 78 and 79 reduce to 


N = A 22 e 0 


and 
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M = - D 
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(k - k q ) (88), (89) 


Boundary Conditions 

Care must be exercised in the application of the boundary condi- 
tions. All boundary conditions must be stipulated in the context of the 
laminate reference surface, as opposed to the laminate midsurface. 

For pinned supports at the middle surface, there is no radial or 
tangential displacement, that is, 


5^ 

1+ 

a> 

o 

V* 

I 

CL 

tt 

o 

(90) 

u(±e Q ,-d) = o 

(91) 


In terms of the reference surface displacement, these conditions are, 
using Eqs. 58 and 59, 

w(±e o ,-d) = w(±e 0 ) = 0 (92) 

u(±e 0 ,~d) = u(±e 0 ) - dft(±e Q ) = o (93) 

Using the definition of ft, Eq. 32, and defining 

= <5 , (94) 

the boundary condition for zero tangential displacement at the middle 
surface becomes 

(1 - 6)u(±e o ) - 5w‘ (+e Q ) = 0 • 


(95) 
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At the pinned supports the bending moment for the middle surface 
vanishes. If V, N”, and M are the force variables for the middle sur- 
faces and V, N, and M are the corresponding force variables for the 
reference surface, the following static equivalence holds. 



V = V V = V (96a) 

N = N N = N (96b) 

M=M+dN M=M-dN (96c) 

The moment-free boundary condition, at the middle surface, is then 

M (±e 0 ) = 0 . (97) 

Referred to the reference surface this condition becomes 

M(±e 0 ) - dN(±e 0 ) = o . 


(98) 
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Using the constitutive relations, Eqs. 88 and 89, Eq. 98 becomes 


D 22^ +e o) ~ K o^ " dA 22 e 0 ^ 0 o^ ” 0 * 


( 99 ) 


Using the definition of k and Eqs. 23 and 41, and defining 


\ a 


J 22 


( 100 ) 


the moment-free boundary condition can be written as 


w"(+.e.) + (l + \)u*(±t\j + w(te_) = o . 


( 101 ) 


By virtue of the vanishing radial displacements at the supports, Eq. 92, 
Eq. 101 simplifies to 

w"(.n* 0 ) + (1 + \)u'(t0 o ) = 0. (102) 


Equations 90, 95, and 102 become the pinned end boundary conditions on 
the kinematic variables u(q) and w(£). 

Transjti on_ Conditions at th e Point Loa d 

Consider the distributed force per unit reference surface area q(9) 
to degenerate to a line load per unit generator direction located at 
0 = ij> as shown in Fig. 5. Define 

Rq ( 0 ) = Pd(d - if’) , (103) 

where «S in this equation is the Dirac delta function. The equilibrium 
equations, Eqs. 60-62, become 



ti If a f* . i> J 
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dV 

de 

dN 

de 

dM 


- N = P<5(9 - tp) cos a 


+ V = P6(6 - sin a 


( 104 ) 


(105) 


+ RV = 0 . (106) 

Integrating these equations between - e and i|> + e, e>o, results in 


ip + e 


dV 

de 


de - 


Ip - £ 


ip + e 


ip + e 
u - e 

Ip + £ 


ip + e 


N d9 = P cos a 


f , 

w de 

A A 

V de = P sin a 

i 

i p - e i 

!> - £ Ip 


6(0 - Ip) de (107) 

Ip - £ 

Ip + £ 

6(6 - ip ) de 

III - £ 


(108) 


ip + e 

r dM 

de 


ip + e 


de + 


RV de = 0 


(109) 


W - £ 


yj - e 

The limit of these equations as e + 0 becomes 
V(i|/ + ) - V(ip~ ) = P cos a 
N(tp + ) - N(ip“ ) = P sin a 
M(/) - M(ijf) = 0 . 


(HO) 

( 111 ) 

(112) 


These same relations can be obtained by considering a free-body diagram 
of a panel segment at the location of the line load. From these relations 


it is clear that the shear and tangential force resultants, V and N, 
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tjf. 


experience a jump at the location of the line load. The moment, however, 
is continuous. The continuity condition on M is similiar to the continuity 
condition on several kinematic quantities which remain continuous at the 
line load. The kinematic quantities which remain continuous are the 
tangential displacement, the radial displacement, and the cross-sectional 
rotation. Specifically, 


w(> ip ) = w(ijT) 


u(4» ) “ u(ljf) 


Q(ip ) = fi(^~) . 


(113) 

(114) 

(115) 


All transition conditions can be stipulated in terms of the funda- 
mental kinematic variables u(e) and w(9) by substitution of previously 
defined quantitites. Working first with those quantities which are 
continuous at the point load, continuity of rotation, Eq. 115, becomes 


w 1 (i|/~) + u(/) _ w 1 (y~) + u(nT) 

R R 


(116) 


Equation 32 was used for this equality. By virtue of the continuity of 
the tangential displacements, Eq. 114, continuity of rotations becomes 

wV) = w'(tff) . (117) 

The jump in the membrane force, Eq. Ill, can be represented with 
the aid of Eq. 88 as 


A 22 {e e (i|! } - = P sin a . 


(118) 
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From the defi niton of e Q , Eq. 23, this equation becomes, with a slight 
rearrangement of terms and using continuity of radial displacements, Eq. 
113, 

u'(f) + ^sin-2L * u '(/) . (119) 

M 22 

By using the constitutive relation between moment and curvature, Eq. 89, 
continuity of moment, Eq. 112, becomes 

-D 22 (k(i|i + ) - k q ) = -D 22 (k(^“) - k q ) . (120) 

From the definition of k, Eq. 41, continuity of moment reduces to, 

wV) + u'(^ + ) = w"(<jT) + u'(^) . (121) 

Finally, the jump in the shear force can be examined. Since shear 
deformations are ignored the shear force must be written in terms of the 
other force variables by using the equilibrium equations. Using Eq. 62, 
the jump in the shear force, Eq. 110, can be stated as 

- I M' (i|» + ) + 1 M' (if) = P cos a . (122) 

Using steps similiar to the previous steps, this condition can be 
written in terms of kinematic quantities as 

3 

+ u"(<jT) + ^n - 9 -— = w "'(i|> + ) + u"(/) . (123) 

u 22 

Equations 113, 114, 117, 119, 121, and 123 are the conditons the kine- 
matic variables u(e) and w(0) must satisfy at the point of application 
of the load. 



32 

04 omesssmu; : r : 

OF POOR 

Summary of Equations 

For convenience the important equations are summarized here. The 
equations, as a whole, constitute a properly formulated boundary value 
problem which can be used to study the linear elastic response of an 
arbitrarily deep cylindrical laminated composite panel to a line load. 

A symmetrically laminated panel is but a special case. 

The equilibrium equations in the intervals -0 <8«{i and ^*<e<0 are: 


V' - N - 0 


(q = 0 for line load) 


(60) 


N' + V = 0 


M' + RV * 0 . 


(q = 0 for line load) 


(61) 

(62) 


The constitutive behavior of the panels is: 


N = A 22 e 0 = A 22 ^ R w 


M = -D 22 (k - k q ) = -D 22 


w" + 


R' 


4 


(88), (23) 

(89), (41) 


The boundary conditions for the pinned supports are: 


radial disp. = 0: 


w(±e Q ) = 0 


(92) 


tangential disp = 0: (1 - 6)u(±0 Q ) - <$w'(±0 o ) = 0 (95) 


moment = 0: 


w"(±6 0 ) + (1 + X)u l (±0 Q ) = 0 


( 102 ) 


At the line load the following transition conditions must be enforced: 
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continuity of radial disp.: 
continuity of tangential disp.: 
continuity of rotations: 
jump in membrane force: u'(\jf) 

continuity of moment: w"(i}T 


w(*“) a w(4) + ) 

013) 

UUO «* u(* + ) 

(114) 

w* (iff) G w' (l|j + ) 

(117) 

^JlSLAn „V) 

(119) 

m 22 


+ u 1 (if") = w"((Ji + ) + u' 

(«' + ) 021) 


jump in shear force: 


w' 


(*•) + u>t*') + B’gSSS - a . (*♦) + „"(/) 


22 


(123) 


By virtue of the nature of the problem, there will be two solutions 
to the homogeneous governing equations given by Eqs. 60-62. One solution 
will represent the radial and tangential displacements of that portion 
of the panel to the left of the point load. The other solution will 
represent the displacements of the panel to the right of the point load. 
These two solutions will be related through the transition conditions 
and each solution will satisfy the boundary conditions pertinent to that 
portion of the panel. 


3. SOLUTION TO GOVERNING EQUATIONS 


General Solution to the Differential Equations 


From Eq. 61 with q = 0 

V = -N' . (124) 

Substituting this into Eqs. 60 and 62 yields 

N" + N = 0 (125) 

M' - RN' = 0 (126) 

By Eqs. 88 and 89 the equations become 

eg + e Q = 0 (127) 

and 

° 22 (k - k 0 )' + RA 22 e 0 ' = 0 (128) 


The solution to Eq. 127 is 

e 0 = A-] cos 0 + A 2 sin 0 . (129) 

Integration of Eq. 128 produces 

D 22 (< - < 0 ) + RA 22 e 0 = A 3 R ‘ 


Using the definition of the reference surface strain, e fl , and curvature, 

u 

k, Eqs. 23 and 41, Eqs. 129 and 130 become 


031) 


u' = w + RA-| cos 0 + RA 2 sin 0 


w" + u' = 


R ^22 / \ 

“ D^“( A 1 cos 6 + A 2 sin 6 ) 


+ X A 

D 22 ”3 ■ 


(132) 


Subtracting the two equations and defining 


*1 


22 
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<Po - R + 


r 3 a 


22 


J 22 


034) 


results in an equation for w(0), namely 


w 


" + w - - A-|4> 2 cos 0 - A 2 4> 2 sin 0 + A 3 


035) 


The solution to trils equation is 

<pp <pp 

w(0) = ~ ~2 ® sin 9 + A 2 0 cos 0 + A^^i + 


A^ cos 0 + A 5 sin 0 


0 36) 


Substituting into Eq. 131 and integrating once leads to the expression 
for u ( 8 ) , specifically, 


u(e) = A 


1 


(«- 


<t>' 


sin e + -s- e cos e 


-(*- 4 ) 


*2 \ . *2 


-£j cos e + 0 sin e 


A^^i 0 


A 4 sin 0 - A 5 cos 0 + Ag . 


(137) 


3S 


: . > 

iff 


Note, Eqs. '36 and 137 are valid in regions of the panel both to the 
left and to the right of the point of application of the line load. 
However, constants thrash Ag will be different for each region. In 
what follows the superscrip of +’s or -'s will be used to denote 
whether the various quantitites being discussed pertain to regions of 
the panel to the right or to the left, respectively, of the line load. 
Thus the solution for w(e) and u(0) to the left of the point load will 
involve constants A^ through Ag while the solution for w(6) and u(e) to 
the right of the point load will involve constants A-| through Ag. 


Application of the Boundary and Transition Conditions 

Since the boundary conditions and transition conditions have been 
written in terms of the kinematic quantities u ( 8 ) and w(e), application 

•f + 

of these conditions to find the twelve unknown constants A, through A g 
and A^ through Ag is straight-forward. In terms of the twelve constants, 
the conditions lead to the following twelve equations: 
radial displacement at 0 = -0 is zero: 


♦l A 3 " | A 5 + “TT A 1 | s1n 0 o + 
radial displacement at 0 = +0 Q is zero: 

* 1 A 3 + | fl 5 - T 2 - A 1 I s1n S o + 




+ 


*2 e o 

2 


Al 


cos 0 Q = 0 

(138) 


COS 0 Q = 


0 

( 139 ) 


37 


o::rn% 


tangential displacement at 3 ® -3 Is zero: 


[HJP% (I 


(1 - 1 ) [a- - - !a: ♦ /r - a: . fifo 


*1\ A “ * 2 a o A -t 
2 A i - — r a e i 


A 1 "TV) sin 0 O 


a: + U - ^ a; * ^ a;'! 


2 - -r A, ( COS e 0 - -s | j a 4 + -f a- 


2 A 2 l s n e o + i A 5 + T A ? + a 'm cos a] » 0 


( 5 2 h 2 


tangential displacement at 8 = +e o is zero: 

(i-s)[a** w ; + |a; + (r.^a;^a 2 + 


j s1n e o 


A 5 + ( R - t) A 2 “ "V A li cos d 0 - 4- *}aJ + % A| 


^ 2 ^o . + ) , ( + *9 . ^, 3 . , ) 

~ A« | sin 9 + j A c + — ~ A i - —^— 2 . a + ’ 

*• 2 ) o | 5 2 *2 2 rt l l 


A-j cos 0 Q = 0 


moment at -3 = -< 3 q i S zero: 


j S + >2 A 2 + “¥ A i j sin e Q - j A" + * 2 A’ - ^ a;[ cos 0 £ 


+ (i + \) - | Ag + ra; + 


“T~ A i I S1n $ 


I - - i>-) 0 n ) 

+ | A 4 + RA 1 " "V" A 2 cos e 0 = 0 


( 142 ) 


< Q 
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moment at 0 * +9„ is zero: 
o 


| A 5 * *2 A 2 - A 1 J sin 3 o - j A 4 * *2 A I + A 2) 


COS 9, 


+ (1 + x) 


♦l A 3 + j A 5 + RA 2 ■ "V' A l'| si " 0 ! 


+ j» + + RA* + AjJ COS e 0 


continuity of radial displacements at 8 ■ ip: 

( . <M - ) ( <M ) 

*1 A 3 + | A 5 - — a t j sin ip + j A“ + A" j cos <p = 

“h A 3 + | A 5 " ~r A lj sin ^ + | A 4 + "r A 2} cos * 

continuity of tangential displacements at 0 = ip: 


(143) 


(144) 




<p~ip 


A 6 +< V A 3 + j A- + ( R F ) A] + ~Y~ A 2 


sin ip 


N 


a: + r - 


V 


2 ) ^2 * 2 j cos tjj = 




jA 5 + /r - r4 )a; - 


^2 \ + ^ 2 ^ +) 

1 Ag * 2 A i | cos ip 


(145) 
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continuity of rotations of 9 = i p: 



Ag sin tp + cos ip + - p — — = a^ sin to + a| cos $ (147) 

continuity of moment at 0 = ip: 

4>i + (R - ^ 2 )^2 sin ^ + (R " 4>2 ) A-j cos ip = 

^1^3 + (R “ 4> 2 )A 2 S ^ n 4* + (R - 4> 2 )a| cos ^ (148) 

jump in shear force at 9 = \p: 

- (R - $2^] sin V + (R " ^2^2 cos ^ + P cos a = 

- (R - 4> 2 )a| sin ip + (R - 4> 2 ) A 2 cos ^ • (149) 

Equations 138-149 are simply a set of linear algebraic equations of the 


form 
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(150) 


^ j ^ J 

These can easily be solved for A^ through Ag. The Appendix gives the 
coefficients C.. and B., i , j = 1,12. 

i J J 

Given the cylinder geometry, the material properties, and the 
loading condition, the matrix C and the vector B can be defined and 

-f. 

coefficients A^ through Ag can be computed. Once these coefficients are 
known, the deformed shape, the strains, the membrane force, the banding 
moment, the boundary reactions, and any other quantity related to the 
response can be determined. The next section presents numerical results 
regarding the resDonse. 
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4. NUMERICAL RESULTS 

Figure 6 shows the predicted stiffness characteristics of three 4- 
layer cylindrical panels as a function of panel geometry. The panels 
are assumed to be fabricated from glass-epoxy with the following mate- 
rial properties: 

E 1 = 20.0 x 10 6 

E 2 = 1.04 x 10 6 

6 12 » 1.50 x 10 6 

v-j 2 = 0.28 

t = 0.005 in. (lamina thickness) 

The case shown in Fig. 6 is for panels with the line load exactly at 
mid-span (ip = 0) and aligned radially (a = 0). The stiffness * K* is 
defined as the total load divided by the radial deflection under the 
load. The stiffness has been normalized by the stiffness of a pinned- 
end panel with a symmetric [90/0/0/90]^ stacking arrangement and with a 
semi -opening angle of tt (see Eq. B.9 in Appendix B). The pinned-end 
panel with semi -opening angle of tt is not be be confused with a closed 
ring. The difference is illustrated below. 
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pinned-end panel with 
semi -opening angle tt 
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closing ring 


The notation on stacking sequence needs some explanation. The 
notation reads from the left to the right, for example in [90/0/0/90]^, 
and each number is associated with a particular lamina. The number 
indicates the angle of the fibers in the lamina relative to the global 
x-axis. The left most entry in the stacking sequence notation refers to 
the lamina at the negative most z position. Reading to the right in the 
notation corresponds to moving from lamina to lamina in the direction of 
positive z. Thus, for the problem here, [90/0/0/90] j signifies that the 
fibers in the innermost and the outermost lamina are in the circumferen- 
tial direction while the fibers in the two inner lamina are in the 
generator direction. 

In Fig. 6 the abscissa is the scaled semi-opening angle 9 0 /tt. 

This parameter completely characterizes the panel's depth, or "bulge," 
with respect to the chord line connecting the support points. It can be 
seen from Fig. 1. that this rise to span ratio for the panel is 



Fig. 6 Stiffness Characteristics of Cylindrical Panels of Various Stacking Arrangements 
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Thus, as stated, the semi-opening angle completely characterizes the 
panel's depth. Moreover, Eq. 151 implies that for a shallow panel in 
which the rise to span is small that 0 Q is small, i.e., 

~ j 9 0> for e Q small . (152) 

Shallowness means a small semi -opening angle, independent of the panel's 
radius. 

Figure 6 shows the stiffness characteristics of a [90/0/0/ 90] T 
panel, a [0/0/90/90] T panel and a [90/90/0/0] T panel. The first panel 
represents a symmetric stacking arrangement that can be made using 
ordinary fabrication techniques, i.e. laying the prepreg over the cylin- 
drical mandrel and curing it. The second panel represents an unsym- 
metric lay-up that would naturally assume the cylindrical shape if it 
were fabricated flat, cured, and cooled. The third panel is a fictious 
unsymmetric panel which if actually fabricated would have curvature 
characteristics opposite to that implied in Fig. 1. This panel was 
studied only to assess the effects of the through-the-thickness location 
of the circumferential fibers. It is important to note that all of the 
panels in Fig. 6, and in subsequent discussions, have an equal volume of 
material (arc-length of 20 inches). Thus, in a sense, material efficiency 
is being studied. 
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From Fig. 6 it appears that the symmetric panel affords the stiff- 
est configuration, except when the opening angle is small. The stiff- 
nesses of the two unsymmetric panels are identical. The small opening 
angle region will be referred to as the shallow panel case. Before this 
study began, it was not clear whether the unsymmetrical ly laminated 
panel would be softer than the symmetric counterpart. If membrane 
forces governed the stiffness, perhaps the stacking arrangement was not 
important. If bending forces governed the stiffness, then the stacking 
arrangement would be important. With a [90/0/0/90] T configuration, the 
circumferential fibers are in the best position, through-the-thickness, 
for resisting bending. In the other stacking arrangements the circum- 
ferential fibers are all on one side or the other of the bending axis. 
Thus it appears bending stiffness is quite important, at least for deep 
panel configurations. 

To obtain a more detailed comparison of the stiffness when shallow 
panels are considered, the stiffness and geometric perimeters are 
rescaled. The stiffness is nondimensional ized by the stiffness, K^, 
of a [90/0/0/90]y simply supported beam (see Eq. B.20 Appendix B) and 
the ratio of the panel rise to laminate thickness is used to describe 
the deviation from a flat configuration. Figure 7 shows a comparison of 
the 3 panels for the case of shallow panels. The figure has several 
interesting features. First, it is obvious a slight curvature in a 
'beam' dramatically increases the stiffness. For a rise-to-thickness 
ratio of even 1, the stiffness has increased several times. This 
phenomena is not unique to composites. It happens with metals and the 
increased stiffness is due to geometric effects. The second feature to 



Fig. 7 Stiffness Characteristics of Shallow Panels 
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be noted is that for a small range of rise-to-thickness ratios, the 
unsymmetric [90/90/0/0]^ panel is stiffer than the symmetric panel. 

Above a rise-to-thickness ratio of about 1, the [90/90/0/0J T is always 
less stiff than the symmetric panel. The ratio of the stiffness for 
each unsymmetric flat panel to the stiffness of the symmetric flat panel 
is 0.545, which may be determined from Eqs. B . 1 9 and B.20 in Appendix B. 

For a shallow panel, the induced membrane force is a significant 
component of the panel's stiffness. This membrane effect adds to the 
inherent bending stiffness. Thus the membrane force at the pinned 
support was examined. For the loading here the membrane force was 
compressive, or a thrust, and Fig. 8 shows how it varies with panel 
geometry for the 3 different stacking arrangements. From the figure it 
is clear that even though the binding stiffness of an unsymmetric panel 
is less than the bending stiffness of the symmetric ones, the induced 
thrust in an unsymmetric panel can be greater. In the case of the 
[90/90/0/0 ]y panel this increased thrust increases the global stiffness 
beyond the symmetric panel's. The figure also indicates that the deeper 
the panel, the less is the difference in thrust among laminates. As 
suggested by Fig. 8, the thrust has significantly less contribution to 
the global stiffness for the deeper panels. 

It is interesting to determine if the shallowness approximation (e 
small) is influenced by the laminate configuration. For metallic panels 
shallowness is determined solely by geometry, see Eg. 152. The error in 
the shallowness approximation on panel stiffness is shown in Fig. 9 as a 
function of the scaled semi-opening angle 6 /t. Panel stiffness as 
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Thrust Versus Rise for Shallow Panels 
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Error in the Shallowness Approximation for Panel Stiffness 
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computed from the assumption of shallowness is designated and is 
derived in Appendix B (see Eqs. B. 1 0 to B.17). The percent error on the 
ordinate of Fig. 9 is defined as 

e = (K sh /K-l)100, 

in which K denotes the exact panel stiffness from the general theory (9 
arbitrary). For the resolution on the graph shown in Fig. 9, the dif- 
ferences between the error curves for the [90/0/0/90] T , [90/90/0/0]^, 
and [0/0/90/90 3 j laminates are indistinguishable. Hence shallowness is 
essentially unaffected by lamination dissymmetry, and thus, it truely is 
a geometric quantity. 

Of fundamental importance to the integrity of composite structures 
is the strain levels. Due to the brittleness of epoxies, large tensile 
strains in the matrix direction of the epoxy matrix laminates cannot be 
tolerated. Even if adjacent layers have the ability to hold the struc- 
ture together, cracks in lamina are simply bad. A crack in one lamina 
can induce cracking in an adjacent lamina even though the cracks could 
ultimately be at right angles to each other. Cracking can lead to 
delamination and the problems inherent with that sort of damage. Thus 
it is important to assess the strain levels in these three laminates. 
Figure 10 shows the circumferential strains in the three laminates. 

The maximum circumferential strains at the innermost and outermost 
radial location were computed as a function of panel geometry. The 
outermost strain's were compressive while the innermost strains were 
tensile. The maximum strain's occurred at the load location. In the 



Fiq. 10 Maximum Circumferential Strain in Cylindrical Panels 
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symmetric laminate, the circumferential fibers in the inner and outer 
plys react these maximum strains. In the [90/90/0/0]^ the inner tensile 
strains are reacted by matrix while the outer compressive strains are 
reacted by fibers. In the [0/0/90/90]^ the tensile strains are reacted 
by the fibers while the compressive strains are reacted by the matrix. 

As can be seen in Fig. 10, the [90/S0/0/0]y laminate yields the 
largest tensile strains while the L0/0/90/90] T laminate produces the 
largest compressive strains. Shallowness appears to have much co do 
wich the magnitude of the strains on the opposite sides of the panel. 

As the panels become deep (H/h increasing) the strains on opposite sides 
become nearly equal in magnitude, an indication of bending effects 
dominating. 

The results presented so far have been for the case where the load 
is actually radial and it is located exactly at midspan. This is an 
ideal situation. In practice it is difficult to achieve the symmetry in 
the loading implied by this situation. Oftentimes the load is slightly 
rotated (a t 0) relative to the radial lines so that it does not act 
entirely in a radial direction. In addition, the load may not be 
exactly at midspan (t{> f 0), acting a little to one side or the other of 
the geometric center of the panel. In some structures such deviations 
from the ideal loading can cause a response that is drastically dif- 
ferent than the response to the ideal case, even if the deviation is 
small. This is particularly true for imperfection sensitive structures 
which exhibit a significant decrease in the static buckling load of the 
perfect structure for small eccentricities in the load position or for 
small geometric imperfections in the structure's initial shape. For 
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example, an elastic arch with symmetric end conditions and unloaded 
shape, subjected to a downward point load at midspan, is imperfect 1 ' on 
sensitive if the arch rise is sufficiently large to cause instability at 
a bifurcation point L5]. Thus the sensitivity of the cylindrical panels, 
particularly the unsymmetrical ly laminated ones, to deviations from the 
ideal loading case was studied. The stiffness of the panels was again 
used as an indicator of panel response. If the load is rotated slightly 
relative to the radial direction, not ar of the deflection is radial. 
However, the stiffness was still defined as the total load divided by 
the radial deflection under the load. The reason for this definition of 
stiffness is that if experimentally one was determining the stiffness 
and did not realize that the load was rotated, this definition would 
still be used to compute stiffness. 

In what follows, the stiffnesses for the nonideal cases have been 
normalized by the stiffness for the ideal [90/0/0/90]y case. The 
stiffnesses will be assumed to be a function of the location of the 
load, x i>, and the inclination of the load, a. Thus, the notation 
< = <(tp,a). The ideal case will be donated as <(ip=0,a=0). Figure 11 
shows the stiffness of a typical shallow panel (R = 100 in., 9 0 /ir = 

0.032) as a function of the location of the load along the span of the 
panel. The load is assumed to act radially (a = 0). Figures depicting 
the effects of load location and orientation can become cluttered if, in 
addition, various panel geometries are shown. Thus a single panel 
geometry representative of a shallow panel was chosen. As can be seen 
in the figure, slight deviations of load locations from the center 
positions cause no noticable decrease in stiffness. As the load location 
reaches quarter-span, 'ii/ Q q = 0.5, the stiffness decreases considerably. 



Fig. 11 Sensitivity of Panel Stiffness to Location of Radial Load (R=l 00 in., 0 /n=0.032) 
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As the point of application of the load approaches the support, iJj/9 q = 

1, the stiffness becomes very large. The stiffness variations of both 
unsymmetric laminates and the symmetric one are similar and there are no 
adverse sensitivity problems with the unsymmetric laminates. In essence, 
if the load were not exactly at midspan in an experiment, the stiffness 
calculations would not be in error. Figure 12 shows the stiffness 
characteristics of a deep panel (R = 6.5 in., 9 /it = 0.5) as a function 
of load location. Again, there does not appear to be any adverse sensi- 
tivity for any laminate. 

Figures 13 and 14 show the effect on stiffness of a rotation of the 
load from the purely radial inclination. Figure 13 is for a represen- 
tative shallow panel (R = 100 in., 9 Q /ir = 0.032) and the load is assumed 
to act at midspan (^ = 0). It is clear that for inclinations 15° - 20° 
from the radial direction, the increase in stiffness over the ideal case 
is negligible. As the load rotates more and more from the radial position 
it begins to encounter the membrane effect of the panel and the radial 
deflections decrease. The unsymmetrically laminated panels seem to have 
the same characteristic as the symmetric panel. In addition, as shown 
in Fig. 14, the deeper panels behave similarly. 



Fig. 12 Sensitivity of Panel Stiffness to Location of Radial Load (R=6.5 in., 0 o /ttO.5) 




Fig. 13 Sensitivity of Panel Stiffness to Load Inclination (R=100 in., 0 /t?=0. 




5. A LIMITED COMPARISON WITH EXPERIMENTAL RESULTS 


To determine if the predictions of the response of unsymmetrically 
laminated cylindrical panels were correct, a [0^/90^]^ T300/5208 graphite- 
epoxy cylindrical panel was fabricated. In the uncured state the flat 
panel was 12 x 12 in. When cured, the panel had a radius of curvature 
of 10.2 in. The panel was fabricated at the Air Force Materials Laboratory, 
Wright-Patterson Air Force Base. The authors extend their appreciation 
to Stephen W. Tsai for providing for the panel fabrication. 

The panel was mounted horizontally, the panel midspan being higher 
than the supported ends. The pin supports were simulated by using two 
steel right-angle members parallel to the generator direction and running 
the lenqth of the cylinder. The open side of each angle member was 
towards the panel. One leg of each angle member was horizontal and the 
other leg was vertical. Each straight edge of the panel was seated in 
the corner of a right-angle. The angles were clamped so they could not 
move relative to one another. The vertical leg of each angle prevented 
the cylindrical panel from spreading as a downward load was applied. 

The horizontal legs provided a vertical reaction to the applied downward 
load. Dead weights acting on a stiff bar that ran the length of the 
cylinder provided the line load. The load acted radially (ct = 0) at 
the panel midspan (9 - 0). The results of Figs. 11-15 indicated that 
the panel response would not be overly sensitive to slightly off-center 
(9 f 0) or slightly misaligned loads (a f 0). Figure 15 shows the 
setup. Dial gages were used to monitor the radial deflections of the 
cylinder at several locations. 
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Figure 16 shows the load deflection characteristics of the panel. 

The experiment was terminated after the panel deflected more than 4 
laminate thicknesses simply because this was felt to be the range of the 
linear theory developed here. The correlation is not particularly good. 
To put the results in context, a [(SO/O^s T300/5208 graphite-epoxy 
cylindrical panel with a 12-inch radius was tested. The panel was 
obtained from the Air Force Flight Dynamics Laboratory and the authors 
extend their appreciation to Capt. Marvin Becker for providing the 
panel. The panel had been previously buckled with inplane compressive 
loads acting in the generator direction. There was no noticeable damage 
to the panel from the previous buckling. Figure 17 shows the symmetric 
panel and the unsymmetric panel side-by-side. Figure 18 shows the load- 
deflection characteristics of the symmetric panel. Again the correlation 
is not good. In fact, the correlation for the symmetric panel is ;.<orse 
than the correlation for the unsymmetric panel. It is not entirely 
clear why the correlation for both panels is not better. Some overstiff- 
ness is expected in the predictions but not to the degree observed here. 



Radial Deflect 
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APPENDIX A 
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Definition of Coefficients C. . and B. 
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cos 9 0 
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C l,2 = ~2“ cos 0 O 
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C 1 ,4 = cos 9 0 
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°1 ,5 

C 1 , 1 2 = 0 


sin 9, 


Eg. 139 


B-j =0 


C 2J ^ C 2 ,6 = 0 

C 2 ,7 = C l,l 
C 2,8 = ‘ C l ,2 
C 2 ,9 = C 1 ,3 
C 2 ,10 = C 1 ,4 

c 2,n = ‘ c i ,5 
C 2,12 = 0 


b 2 = 0 
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Panel Stiffness by Castigliano's Theorem 

The global panel stiffness may be obtained with relatively few 
algebraic manipulations by Castigliano's second theorem for some special 
geometric configurations if the line load P is directed radially inward 
at midspan (ip = a = 0 in Fig. 4). These manipulations are outlined in 
this appendix for the deepest panel (0 = it), a very shallow panel (0 O 

smal 1 ) , and a flat panel . 

The complementary strain energy per unit length in the x-direction 
of the panel is 


V = 


< 7 > 


-0. 


J 22 


22 


]R d6 


(B.l) 


From appropriate free body diagrams of a circular cylindrical panel with 
arbitrary 0 , the statically admissible force fields as a function of 0 
are determined. These are 


N(0) = - T cos 9 - (j^) sgn (9) sin 9 (B.2) 

M(0) = - T[R cos 0 - a cos (9 )] - (j)[R sgn (9) sin - a sin 9 ], 

(B. 3) 

in which |9| < 0 , T is the horizontal reaction thrust per unit length 
in the x-direction exerted on the panel by the pinned-end supports and 


sgn(9) = { 


(B.4) 


1 , i f 9 > 0 
-1, if 9 < 0 . 
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For the deepest panel 0 Q = it. Using this specific value of 0 Q , 
Eqs. (B.2) and (B.3) are substituted into Eq. (B.l) and integrations on 
0 are performed. The resulting complementary strain energy is 


v = Ml 


2A 


(Rtt/ 4)P 2 + R 3 [1 + 2(a/R) 2 lT 2 + (R 3 tt/ 4)P 2 + (4aR 2 )TP 


22 


2D, 


22 


(B.5) 


Castigl i ano ' s second theorem states that partial derivatives of V with 
respect to the forces T and P are the generalized displacements. In 
connection with the specific problem here, this requires 


3 V „ 
ap “ 


A 


and 


H- = 0 , (B.6),(B.7) 


where A is the panel displacement in the direction of the load at the 
point of load application. The partial derivative with respect to T 
(B.7) must vanish since the pinned-end supports are immovable. Eqs. (B.6) 
and (B.7) result in linear algebraic equations for T and P when 
Eq. (B.5) is used for V. If T is eliminated between them, a single 
equation of the form P = K^A results, in which is the global panel 
stiffness. After some manipulations this stiffness is 

K = (4 1 rD 22 /R 3 )[l + 2(a/R) 2 + D 22 /A 22 R 2 )] 

• {tt 2 [1 + D 22 /(A 22 R 2 ][1 + D 22 /(A 22 R ) + 2(a/r) 2 ] - 16(a/R) 2 } 1 . 

(B.8) 

Specializing this for a symmetric laminate (d = 0, D 22 = D^, etc.), and 

p 

noting that the factor D 22 /(A 22 R ) is small with respect to unity and may 
be neglected for the panels in this study, the stiffness in Eq. (B.8) 
is adequately approximated by 
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(B.9) 


* 12uD 29 /[a 3 (3Tr 2 - 16)] . 

This is the normalization constant for the ordinate in Fig. 6. It is 
worth stating that neglecting the factor D 22 / (^22 R ) res P ect t0 one 
is equivalent to neglecting the membrane portion of the complementary 
strain energy with respect to the bending portion; i.e. let A 22 -*■ 05 in 
Eq. (B. 5) . 

Consider a shallow panel which is quantified by a small semi- 
opening angle 0 Q . Statically admissible force fields for tnis approxi- 
mation are 

N(0) = - T (B. 10) 

M(0) = - (H - R0 2 /2)T - [sqn (0)R - a0 Q ](P/2) , (B.ll) 

in which 

H = d + H . (B. 12) 

p 

The rise H (see Fig. 1) in Eq. (B.12) is closely approximated by 30 q/ 2 
for a shallow panel. Substituting the functional forms of the force 
fields into the complementary strain energy (B.l), performing the inte- 
grations on 0, results in 

V = (R0 q /D 22 [C t T 2 - C Tp TP + C p P 2 ] , (B. 13) 

in which 

C T = H 2 - HR0 2 /3 + R 2 0 4 /2O + 0 2 2 /A 22 (B.14) 

C Tp = Ha0 o - HR0 o /2 - aR 3 /6 + R 2 8 3 /8 


(B.15) 
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C p = (a 2 9 2 + R 2 9 2 /3 - aR0 2 )4 . (B.16) 

Castigliano's second theorem is applied to the complementary strain 
energy in Eq. (B . 1 3 ) , and the same procedure used to obtain the stiff- 
ness for the deepest panel described previously is followed. The 
resulting expression for the shallow panel stiffness is 

K sh = (D 22 /Re 0 )2C T /(4C T C p - C 2 p ) . (B. 17) 

Finally, consider a flat panel or wide beam. The results for the 
shallow panel may be used to obtain the global stiffness for this case 
if the following limiting process is undertaken: Let 9 Q 0, such that 

a0 Q -*• L/2, where L is the beam length (see Fig. 1). In this limit, 
then, the coefficients in Eqs. (B.14) to (B.16) become 

C T = D + D 22 /A 22 , 

C Tp = dL/4 , and (B. 18) 

C p = L 2 / 48 . 

Substituting these limiting values into the global stiffness expression 
(B.17) results in a formula for the flat panel stiffness which is 
designated . After some manipulation this stiffness is 

K bu = ( 48D 22 /l3) (d2 + D 22 /A 22 )/((d/2) 2 + D 22 /A 22 ) . (B.19) 

For a symmetric flat panel this becomes 

K b = 458 22 /l * 3 ’ (B,2 °) 


which is the normalization constant used for the ordinate in Fig. 7. 


